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Abstract: This paper gives a short summary of some of the aspects of charmonium 
■ -■ . which can be addressed at BES and CLEO-c and other e~^e~ facihties. These topics 

Q ! include the spectroscopy of charmonium states, radiative transitions, e+e~ widths, 

O ' two-photon widths, hadron loop effects and open-flavor strong decays. 

1^ ' 1 Introduction 

0\ ', Charmonium has been called the "hydrogen atom of QCD", since many of the most char- 

I acteristic and interesting aspects of QCD can be inferred from studies of the spectrum of 

charmonium states and their decays and interactions. Charmonium is a useful system for 
the study of forces between quarks in QCD, especially the poorly understood, nonperturba- 

^ . tive confining interaction. This system is of special interest because the dynamics are only 

^ ! quasirelativistic, and the quark-gluon coupling at this mass scale is intermediate in strength; 

Q I for this reason the more unusual features of QCD, such as relativistic corrections to the 

^ ' properties of bound states, spin-dependent forces from quark motion, virtual meson decay 

loops and various other effects may be only moderately large "controlled" corrections to a 
simple nonrelativistic potential model picture. Thus in charmonium we have a laboratory in 
, which various novel dynamical effects may simultaneously be large enough to be clearly iden- 

D I tified, and small enough so that they can be treated as perturbations of a familiar quantum 

i-^ ' mechanical model. 

^ ■ In this contribution we will discuss some of the aspects of charmonium that are most 

^ ■ easily accessible at BES and CLEO-c. We will primarily discuss charmonium states above 

the open-charm threshold of 3.73 GeV, since these are quite poorly known; it should be 
straightforward to improve our knowledge in this area considerably. For our discussion we 
use a "generic" nonrelativistic potential model to generate the spectrum of states, and the 
resulting wavefunctions are used to evaluate electromagnetic transitions involving low-lying 
charmonium levels. We then discuss open-flavor strong decays of charmonium states, which 
is an important and largely unexplored topic. In our discussion we consider the properties 
of 1 states in particular, as these are the most easily accessed at e~^e~ facilities. 

2 A Generic Charmonium Potential Model 

There are several charmonium potential models in the literature, which for the most part 
share the common features of short-ranged color Coulomb and long-ranged linear scalar 
confining potentials. In addition to this standard "Coulomb plus linear" potential, we also 
include a Gaussian-smeared spin-spin hyperfine term in our "zeroth-order potential" , since 
a careful treatment of the spin-spin force is important for the determination of the ^Si - ^Sq 
mass splittings. This issue has attracted recent attention due to the discovery of the rjc ' by 
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BELLE [T at a surprisingly high mass of 3654 ± 6 ± 8 MeV (their average including inclusive 
production is 3639 ± 6 MeV), since confirmed by BABAR (at 3630.8 ± 3.4 ± 1.0 MeV) 
and CLEO [2j (at 3642.9 ± 3.1 ± 1.5 MeV). The zeroth-order potential model Hamiltonian 
we assume is of the form 

Ho = -— + 6r+ ^^(r)Se-Sg (1) 

rric or 9m^ 

where ^o-(r) = (ct^/tt'^/^) exp (— a^r^). Solution of the Schrodinger equation with this Hq 
gives our zeroth-order charmonium wavefunctions. Splittings within an orbital multiplet 
are then determined by taking the matrix element of the spin-dependent Hamiltonian Hj 
between these zeroth-order wavefunctions. The spin-dependent Hamiltonian is taken from 
the one-gluon-exchange Breit-Fermi Hamiltonian (which gives spin-orbit and tensor terms) 
and an inverted spin-orbit (Thomas precession) term, which follows from the assumption of 
a Lorentz scalar confining interaction. This Hi is explicitly 

jj — "^ i^ g I "^ ^ L S (2) 

The assumption of a Lorentz scalar confining interaction is of course questionable; the nature 
of the spin-dependent contributions of the confining interaction to the interquark force is one 
of the interesting open questions which we hope to address in the study of charmonium. 

There are four parameters in this simple potential model, the one-gluon-exchange cou- 
pling constant a^ (here for simplicity taken to be a constant), the string tension b, the charm 
quark mass rric, and the contact hyperfine interaction smearing parameter a. We have fixed 
these using an equal- weight fit to the the masses of the 11 reasonably well established exper- 
imental charmonium states; the results are as = 0.5461, b = 0.1425 GeV^, rric = 1.4794 GeV 
and a = 1.0946 GeV, which gives an rms error of 13.6 MeV. (Of course no special meaning 
should be attached to these precise values, since the model is somewhat ad hoc and is at best 
only an approximate description of the physics of charmonium.) The experimental states 
used as input and the predictions for the cc spectrum (states below 4.2 GeV and all known 
1 states) are given in Table |21 these and some additional levels are shown in Fig.l. 

It is evident in Fig.l that the simple Coulomb plus linear potential model gives a very 
satisfactory fit to the known charmonium states. It is remarkable that both the mean level 
positions and the splittings within the IS, IP and 2S multiplets are well described, since 
they are testing different aspects of the interquark forces. The mean multiplet positions 
test the accuracy of the "funnel-shaped" potential, and determine approximate values for 
the two parameters «« and b. The multiplet splittings in contrast test whether the spin- 
dependent interactions (spin-spin, spin-orbit and tensor) are also well described with the 
same parameters. These terms arise from the assumptions of transverse one-gluon exchange 
at short distances (which provides the J/i/j-r)c and ip '-rj^ ' spin-spin splitting and dominates 
the P-wave splitting), and scalar confinement (which provides a negative spin-orbit term 
that is required for good agreement with the experimental xj rnasses). The general good 
agreement between the potential model and experiment is remarkable in view of the addi- 
tional, presumably large effects such as virtual decay loops that have not been included in 
the model. 




FIGURE 1. The spectrum of charnionium states predicted by the cc potential model de- 
scribed in the text, together with the 11 known experimental states. The X(3872) is also 
shown, although its identification with cc appears dubious. 

We note in passing that lattice QCD has also been applied to charnionium, and predicts 
a very similar spectrum. (See for example the recent work of Liao and Manke j3].) It is 
very encouraging that the theoretically better-justified LQCD approach gives a pattern of 
energy levels which is very similar to the potential model; this supports our use of the more 
intuitive potential model to describe transitions involving charnionium states in this paper. 
We shall note however in our discussion of decays that virtual decay loops are expected to 
have important effects on both the spectrum and the composition of charnionium states |5|6j . 
and these effects are neglected in pure cc potential models as well as in quenched LQCD. 

3 Electromagnetic couplings 

3.1 El radiative transitions 

Radiative transitions are a very interesting feature of charnionium physics. They are quite 
straightforward to evaluate in cc potential models, and (with sufficient statistics) provide a 
route from the initial 1 states produced in e~^e~ annihilation to C = (+) charmonia. 
The largest rates are for El (electric dipole) transitions, which we calculate using 
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Here Cc = 2/3 is the c-quark charge in units of |e|, a is the fine-structure constant, E^ is the 
photon energy, and the angular matrix element Cfi is 




C;i = max(L, L')(2J' + 1)| j ^1/ " ^^^ 

This formula is the same as that quoted by Ref . ^ , except for our inclusion of a relativistic 
phase space factor of Ej /M^- ^ . (This is usually not far from unity for the cases we consider.) 
To generate numerical results for transitions of particular interest to BES and CLEO-c, 
we have evaluated the matrix elements {n' ^ Lj,| r |n^^+^Lj) using the nonrelativistic 
Schrodinger wavefunctions of the model described in the previous section. The resulting El 
rates are given in Table |3J together with the current experimental numbers. 

Several very interesting features of El radiative transitions are evident in Tabled First, 
the IP -^ IS transitions are in very reasonable agreement with experiment. It is notable 
that the predicted radiative partial width for he — > •yrjc is especially large, which suggests 
that decays to 7?7c may provide a discovery channel for the elusive h^. One possibility is 
77 -^ Vc', followed by the decay chain rjc' — >■ ■jhc, he — >■ 7^c- Since these electromagnetic 
couplings are all reasonably well understood, detection of the h^ would simply be a matter of 
accumulating adequate statistics at a high-energy e+e~ facihty with sufficient 77 luminosity 
at v^ = 3.7 GeV. 

The theoretical rates for the 2S -^ IP transitions appear too large by about a factor of 
two, although the relativized model of Godfrey and Isgur JH] does not share this difficulty. 
We note in passing that good agreement between a pure-cc charmonium potential model and 
experiment may be spurious; decay loop effects will contribute two-meson continuum com- 
ponents to all these charmonium resonances, which may significantly modify the predicted 
radiative transition rates. 

El radiative transitions from the higher-mass charmonium states are especially inter- 
esting. The l^Di candidate iIj{3770) is predicted to have large partial widths to 7x1 and 
7X0 (with branching fractions of 0.5% and 1.7% respectively), but the branching fraction 
to 7x2 is predicted to be only about 2 ■ 10~^. This small number however follows from the 
assumption that the V'(3770) is a pure ^Di state; if there is a significant admixture of S-wave 
basis states in the -^(3770), 

1^(3770)) = cos(e) I^Di) + sin(e) |23Si) (5) 

one typically finds a much larger radiative width to 7x2 j9|10j - (See Fig. |21 The sign of the 
mixing angle 6 depends on the convention for the normalization of the ^Di and 2^Si basis 
states; note for example in Fig.l of ReffTUj that a zero '?/'(3770) -^ 7x2 width requires a 
small negative mixing angle, whereas with our conventions it would be positive.) Since the 
coupling of the ?/'(3770) to e+e~ suggests a significant 2^Si component, a measurement of 
this radiative partial width will be especially useful as an independent test of the presence 
of this amplitude in the ?/'(3770) wavefunction. 

We note in passing that if the dominant mechanism of ^Di - 2'^Si basis state mixing in 
the ?/'(3686) and ^(3770) is through virtual charm meson decay loops such as DD, DD* and 
D*D*, the assumption of a 2 ® 2 orthogonal mixing matrix as in Eql^is incorrect. In this 
case the (^Di|'?/^(3686)) and (2^Si|'?/^(3770)) overlaps will no longer be simply related, and 
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nominal D^ quark model state 



FIGURE 2. Predicted radiative partial width for the El transition -(/'(STTO) -^ 7x2 as a 
function of the 2^Si-^Di mixing angle 9. 



radiative transition amplitudes will also receive contributions from photon emission from 
the two-meson continua. 

Next we consider radiative decays of the higher-mass vectors ■?/'(4040) and ■?/'(4159). As is 
evident in TablelH if the '?/'(4040) is dominantly a 3^Si state as we assume here, it should have 
very small El radiative widths to the triplet members of the IP multiplet, with branching 
fractions of at most about 10~^. The radiative widths to the unknown 2P triplet states 
are theoretically much larger, with branching fractions of 0.3-1. -10"^. Since the 2P states 
have large strong branching fractions to DD {xo and X2) and DD* {xi and X2) [HJ, it 
may be possible to identify these states in the DD and DD* invariant mass distributions of 
^'(4040) -^ 7DD and 7DD* decays. 

Radiative decays of the ^(4159) share certain features with the decays of both the -0(3770) 
and the V'(4040). As we found for the ?/'(4040), the El coupling of the ^(4159) to the 2P 
multiplet is much stronger than to the IP multiplet, so radiative decays of the iJj{4:159) could 
be used to search for 2P states. The branching fraction to the 2'^Po Xo state in particular 
is predicted to be a relatively large 0.6%. A strong suppression of '?/'(4159) El decays to 
n^P2 states is predicted, as we found previously for the ip{3770). This again follows from 
the assumption that the state is purely D-wave cc. Given an important S-wave component, 
which is suggested by the e+e~ width, the coupling of the '?/'(4159) to n'^P2 states should be 
much larger. Finally, a new feature is that one may reach the F-wave multiplet (specifically 



the ^F2 state) through the radiative decay of the '?/'(4159). Unfortunately this ^F2 state is 
expected to be rather broad; we anticipate a total width of about 160 MeV, mainly to DD 
p. 

3.2 Ml radiative transitions 

Ml transitions between charmonium states in pure cc models result from photon emission 
through the Hi = —jl- B magnetic moment interaction of the c quark (and antiquark), and 
as such are suppressed relative to El transitions by the small factor of l/rUc in the magnetic 
moment operator. The Ml transition amplitude is proportional to the matrix element of 
the spin operator, with a spatial factor that (without recoil corrections) is simply the matrix 
element of unity. Ml transitions are therefore nonzero only between states with the same L^g 
(and different Sec, since the C-parity must change). If we assumed a spin-independent zeroth- 
order potential and neglect recoil effects, Ml transitions between different radial multiplets 
would vanish because the n^Si and n'^So states have orthogonal spatial wavef unctions. One 
such transition is actually observed in charmonium, ip' -^ 'yrjc, which must be due in part 
to the nonorthogonal ip ' and rjc spatial wavefunctions and final meson recoil effects. The 
formula for Ml decay rates analogous to the El formula used in the previous section is 



Em, (n --^L, -. n' -'^^L', + 7) = ^.^ ;^ E^ ^ ^ 5,,, 5s,s'.i 



m2 ^^ M^'^^ 2L + 1 



■|(n'2S'+iL'j,|n2S+iLj)|2. (6) 

Evaluating this formula for transitions from the ip and ip ' gives the results shown in Tabled 
A more detailed study of Ml radiative decay rates, incorporating recoil corrections (which 
are numerically important for transitions between multiplets such as 2S ^ IS) will appear 
in Ref.m]. 

TABLE 1. Theoretical and experimental Ml radiative partial widths of the tp and tp', 
neglecting recoil effects. 



Initial state 


Final state 


Tthy. (keV) 


Texpt. (keV) 


J/^ 


IVc 


2.9 


1.2 ±0.3 


r 


IVc' 


0.21 




r 


7Vc 


4.6 


0.8 ±0.2 


Vc' 


7J/^ 


7.9 





A well-known problem is evident in the decay rate J/ip — > jrjc, which is that the predicted 
rate in the nonrelativistic potential model is about a factor of 2-3 larger than experiment. 
Since this rate only involves the charm quark magnetic moment, and hence only its mass, 
this discrepancy is a surprise. The relativized Godfrey-Isgur model jS] predicts a somewhat 
smaller rate of 2.4 keV, which is still about a factor of two larger than experiment. Since 
the errors are rather large, it would clearly be very interesting to improve the experimental 
accuracy of this surprising partial width. If this discrepancy is confirmed, it may be an 
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indication that pure-cc models are a rather inaccurate description of charmonium, and that 
other components of the state vector such as two- meson continua make comparable important 
contributions to the Ml transition amplitudes. In view of the inaccuracy of the theoretical 
IS Ml transition rate, it would also be interesting to test the 2S transition rate ip ' — >■ 7?7c ' 
experimentally. Unfortunately, this rate is predicted to be a rather small 0.21 keV in the 
nonrelativistic cc model. 

An even larger discrepancy between experiment and theory is evident in the "hindered" 
Ml transition t/^' — »• 777c- Since this rate is only nonzero due to recoil effects (not included 
here) and corrections to the naively orthogonal IS and 2S cc wavef unctions, the discrepancy 
is perhaps less surprising than that found in the allowed IS — ;> IS V' — >■ 777c transition rate. In 
any case this is another example of an Ml transition rate in charmonium in which experiment 
and theory are clearly in disagreement. Since the experimental rate is again only about 4cr 
from zero, it would be very useful to improve the accuracy of this measurement at BES or 
CLEO-c. 

Ml decays between charmonium resonances have only been observed between S-wave 
states. The rates between orbitally-excited states are typically predicted to be quite small, 
due to the small splittings within excited-L multiplets. They are large enough however to be 
observable given narrow initial states and large event samples. For example, a hypothetical 
^D2 cc assignment for the X(3872) could be tested through a search for its Ml decay to 
7'?/'(3770), which has a partial width of 1.2 keV in our nonrelativistic potential model. The 
he (assumed at 3525 MeV) decay h^ -^ 7x0 has similar phase space, and is predicted to have 
a partial width of 0.8 keV. In contrast, the smaller phase space of the Ml transition from the 
higher-mass Xi state leads to an expected partial width for X2 -^ 7^c of only about 60 eV. 

3.3 e'^e" widths of 1 states 

Leptonic partial widths are trivially accessible at e~^e~ machines, and they provide interesting 
(and currently rather puzzling) information regarding the wavefunctions of 1 charmonium 
states. In the nonrelativistic cc potential model e~^e~ annihilation through the photon couples 
to the cc pair at contact, so these widths are only predicted to be nonzero for S-wave cc 
systems. This nonrelativistic partial width is given by the van Royen - Weisskopf formula 

Tec- (Si)-ya ^^ (7) 

where the radial wavefunction is normalized to Jl^ r'^dr\ilj{r)\'^ = 1. Strictly speaking, for 
relativistic bound states the annihilation is not completely local, so some coupling of D-wave 
states to e~^e~ is predicted. This width at leading order is proportional to |'?/'"(0)p. 



^cc ( l-'l) = TT" 



^^^.„_ , 50 2 \rm' 



9 Ml,mt 

which is much smaller than the corresponding widths of the ^Si states. Evaluating these 
widths using the nonrelativistic quark model wavefunctions described here gives the results 
quoted in Table |21 

Note that the agreement between the model and experiment is not especially good. This 
may be in part due to the choice of parameters, although the overestimate of the J/'?/' 
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TABLE 2. Predictions of the nonrelativistic cc potential model for e~^e partial widths, 
together with current PDG experimental values |13j . 



State 


Asst. 


r?:r (keV): 


r^;\ (keV): 


J/V' 


l^^Si 


12.13 


5.40 ±0.17 


^' 


2 3Si 


5.03 


2.12 ±0.12 


^(3770) 


l^Di 


0.056 


0.26 ±0.04 


V^(4040) 


S^Si 


3.48 


0.75 ±0.15 


V^(4159) 


23Di 


0.096 


0.77 ±0.23 


ij{U15) 


43Si 


2.63 


0.47 ±0.10 



leptonic width by roughly a factor of two seems to be a common difficulty in naive potential 
models. Non-valence components in the charmonium states, such as ccg or D meson pairs 
from decay loops, may also contribute to this inaccuracy. Concerns have also been expressed 
that the leading-order QCD radiative corrections may be large ^3], which could significantly 
reduce the overall scale of the leptonic widths. Of course the leading-order corrections are 
prescript ion- dependent, so it is not clear that this claim is reliable. In any case it is evident 
that a simple change of scale will not resolve the discrepancies in Table El since the nominally 
D-wave states '?/'(3770) and ip{4159) both have significant leptonic widths. 

The large experimental leptonic widths of the '?/'(3770) and ■?/'(4159) relative to predic- 
tions for pure D-wave cc states (see Table Ej) may be due to an admixture of an S-wave cc 
component |9|10j . As we noted previously, this mixing can be tested through measurements 
of the El radiative decay rates 'i/'(3770) -^ 7x2 and (with more difficulty) '?/'(4159) — *> 7x2 
and 7x2 ', since these transitions are very sensitive to the presence of S-wave cc components. 
Radiative width ratios such as T{ip{3770) -^ 7X2)/r('?/'(3770) — > 7X1) and the leptonic width 
ratio r(-?/'(3770) -^ e+e^)/r('(/;' -^ e+e^) provide two independent tests of S-D mixing. 

3.4 Two-photon couplings 

Although neither BES nor CLEO-c will have adequate ^/s to usefully exploit two-photon 
production of charmonium resonances, this has been a useful technique at higher-energy 
e~^e~ facilities, which we briefly mention here in the interest of completeness. This subject 
has been reviewed elsewhere [T^ . 

Two-photon resonance production involves the reaction e^e^ — > e'^e~R, R —^ f, where 
i? is a C = (±) meson resonance and / is the exclusive flnal state observed. This process 
proceeds dominantly through the two-photon coupling of the resonance, e~^e~ -^ e"'"e~77, 
77 — >■ R, and hence implicitly gives the two-photon partial width r^^(i?) times the branching 
fraction B{R —>■ /). One may also consider the case of one or both photons signiflcantly off 
mass shell g^ 7^ 0, which gives the generalized widths T^*-y{R) and T^*^*{R). The interesting 
question for off-shell widths is whether they are given by a vector dominance formula with 
the mass of the relevant vector, in this case the J/ip. 

In the limit of large quark mass (and hence a zero-range charm quark propagator) the two- 
photon width of an S-wave charmonium state is proportional to the wavefunction at contact 



squared. In this approximation, higher-L states are produced with amphtudes proportional 
to the L*'* derivative of the wavefunction at contact jTH]. (In practice this attractively simple 
contact approximation appears marginal at best at the cc mass scale.) Both the r]c and rjc' 
have been seen in two-photon collisions; the 77 width of the rjc is approximately 7 keV, 
comparable to quark model expectations. The rjc' — > 77 width is not known because one 
measures the 77 width times the branching fraction to an exclusive final state, and the r/c ' 
absolute branching fractions are not known. Higher-L cc two-photon widths are suppressed 
by the mass of the charm quark, and therefore are not very well established; only the P-wave 
Xo and X2 states have been observed. The theoretical ratio of P-wave 77 widths in the large 
quark mass limit is r^^(^Po)/r^-y(^P2) = 15/4, however this may be modified significantly 
by QCD radiative corrections and finite quark mass effects. Nonrelativistically the 77 state 
produced by the ^P2 should be pure helicity two; relativisitic corrections are expected to give 
rise to a small helicity-zero 77 component as well. 

Although one can in principle produce higher-L cc states in 77 collisions, in practice 
these rates fall rapidly with increasing L for cc and heavier QQ systems. As an example, the 
two-photon width of a hypothetical ^D2(3840) cc state is predicted to be only 20 eV [T3] . 

4 Open-charm strong decays 

The dominant decay modes of most hadrons are open-flavor strong decays, in which the 
initial quarks separate into different final states. Despite their importance in QCD, these 
open-flavor decays are not at all well understood on a fundamental level. There are very 
exciting prospects for greatly improving our understanding of this important aspect of non- 
perturbative QCD through studies of the open-flavor decays of charmonium states above the 
DD threshold. 

Strong decay couplings are also of great interest because at second-order they give rise 
to {cc) -^ (cq) (qc) — > (cc) loop diagrams, which are not included in pure cc potential models 
or in quenched lattice gauge theory. As these virtual decay loop effects are believed to 
be numerically important in charmonium |SE], it is correspondingly important to improve 
our understanding of the strong decay couplings that drive these processes. This will be 
most straightforward through studies of the decay couplings of cc states above open-charm 
threshold. 

At present these open-flavor decays are treated using simple quark model descriptions 
of the decay process, which are known to give numerically reasonable results for many light 
hadron decays. The earliest of these decay models is the "^Pq model" of Micu |I2| and 
LeYaouanc et al. ^H] , who describe the decay in terms of the local production of a new qq 
with 0"''+ quantum numbers (hence ^Pq)- The overall strength of this dimensionless pair 
production amplitude 7 is simply fitted to experiment, and for light mesons is found to be 

7^0.4 [in]. 

The best known theoretical work on charmonium is the series of papers by the Cornell 
group 120121122^3] , which includes a study of open-flavor strong decays |23] ■ Remarkably, 
the strong decay model assumed by the Cornell group is not the standard ^Pq model. Instead, 
they assume that strong decays take place through pair production from the linear confining 
potential, which they assume to transform as the time component of a Lorentz vector, 70. 
(Vector confinement was assumed in some hadron models in the 1980s, but has since been 



largely replaced by the assumption of Lorentz scalar confinement.) 

A more recent study of possible QCD strong decay mechanisms by Ackleh et al. 
found that the ^Pq model decay amplitudes were quite similar in form to what one would 
find for pair production by the linear scalar confining interaction, although the usual string 
tension leads to an overestimated decay amplitude by about a factor of 2-3 for light meson 
decays. 

These two models of the decay interaction, Lorentz vector versus Lorentz scalar, will give 
different predictions for open-flavor strong decay amplitudes that can be tested at BES and 
CLEO-c. Unfortunately the Cornell group did not give detailed amplitude decompositions 
for their Lorentz vector decay model; instead they evaluated a "resolvent" that gave the 
cross section ratio i?, and showed graphical results for R and the branching fractions. The 
required evaluation of decay amplitudes in this model is in progress [23]. Similarly, the 
corresponding ^Pq decay model amplitudes have (remarkably) never been evaluated in detail 
for the higher charmonium states. This work has now largely been completed JJ , and some 
of the results are presented here. In particular we will discuss results speciflc to the easily 
accessible 1 states. 

The results we find for the open-charm decay modes of the four accessible 1 states in 
the ^Pq model, with SHO wavefunctions and parameters of /3 = 0.5 GeV (length scale) and 
7 = 0.4 (pair production amplitude) are shown in Table El The lightest state is the '?/'(3770), 
which is usually considered a ^Di state, can only decay to DD; it is predicted to have a total 
width of 43 MeV, about twice the experimental value. The total widths of the nominally 
3^Si '?/'(4040) and 2^Di '?/'(4159) are both predicted to be rather near the PDG experimental 
values ll3j; 74 MeV (theory) versus 52± 10 MeV (experiment) for the ^'(4040), and 73 MeV 
(theory) versus 78 ± 20 MeV (experiment) for the '?/^(4159). 

The branching fractions of the ?/'(4040) are a famous problem; although phase space 
favors DD over DD*, and D*D* has essentially no phase space, the existing data actually 
suggests that the relative branching fractions are D*D* >>DD* >> DD. (The ratios quoted 
by Mark I at SLAC were 32.0 ± 12.0 : 1 : 0.05 ± 0.03 [HESjO The ^Pg model is at least 
partially in agreement here, as is the Cornell decay model [22]! since both predict a very 
small DD mode. (In the decay models this is a result of a node in the decay amplitude 
of the radially-excited 3^Si state.) The Cornell decay model however predict comparable 
branching fractions to DD* and D*D*. The '^Pq model with our parameters finds a similar 
result (see Table [S]), although there is a node in the DD* decay amplitude that can also 
significantly reduce the DD* branching fraction JT] . It would be very interesting to measure 
and recalculate these branching fractions carefully to determine whether there is indeed 
disagreement with theoretical expectations. Note also that the D*D* mode is especially 
interesting, in that there are three amplitudes in this final state, ^Pi, ^Pi and ^Fi. The ^Pq 
model predicts that decays from S-wave cc states have zero ^Fi amplitudes, and that the 
ratio of the two P-wave amplitudes is ^Pi/^Pi = — 2\/5. 

The 2^Di candidate '?/'(4159) has the same modes accessible as the '?/'(4040), but is pre- 
dicted to have a very different pattern of branching fractions and decay amplitudes. The 
mode DD* is theoretically strongly suppressed, DD should be large, and D*D* should be 
the leading mode. The relative amplitudes within the D*D* final state are predicted to have 
a very characteristic pattern; ^Fi is predicted to be dominant, and ^Pi is predicted to be 
larger than ^Pi, in the ratio ^Pi/^Pi = — v5. 
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Finally, the ?/'(4415) is assigned to 4^Si in potential models. This highly radially-excited 
level has several decay amplitude nodes that make the predicted decay branching fractions 
less reliable. The expected dominant mode among S-wave pairs is D*D*, with the same 
amplitude ratios quoted above for the 3^Si state. 

An interesting new feature of ^(4415) decays is that there is sufficient phase space for 
decays to S+P D-meson pairs; in particular, the interesting new CLEO states D*j(2317) 
and Dsi(2460) can both be made in S-wave (assuming that the Dsj(2317)) is a scalar), in 
^p{AA15) -^ D*D*j(2317) and ?/^(4415) -^ D,Dsi(2460) respectively. (Nominally these decay 
modes are just above threshold, but they are accessible once the width of the '?/'(4415) is 
incorporated.) Since the branching fractions expected to these final states are not especially 
small pi], the '?/'(4415) could provide a straightforward source of a large sample of D*j(2317) 
and Dsi(2460) events. 

5 Conclusions 

In this note we have reviewed some topics of current interest in charmonium which can be 
addressed at the BES and CLEO-c facilities. The topics discussed here were electromagnetic 
couplings and decays (El, Ml, e~^e~ and two-photon), hadron decay loops and open-fiavor 
strong decays. Studies of properties of the ip{3770), ^(4040) and ?/'(4159) are especially 
interesting in this regard. The El decays of the '0(3770), in particular to 7x2, can be 
compared with the e~^e~ width to allow tests of S-D mixing, which may be due to open- 
fiavor decay loops. The relative strong branching fractions of the ■?/'(4040) and the iIj{A159) to 
DD, DD* and D*D*, and in particular the relative size of the three D*D* decay amplitudes, 
allow sensitive tests of strong decay models. Finally, the high-mass tail of the '?/'(4415) may be 
useful as a source of events containing the new narrow resonances D*j(2317) and Dsi(2460). 
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TABLE 3. Predictions of the nonrelativistic cc potential model discussed in the text for the 
charmonium spectrum to 4.2 GeV (and for all S-wave states to 4.5 GeV). 



State Mass (MeV): Expt. (input) 


Theor. 


J/V^ 


3097 


3090 


Vc 


2979 


2982 


r 


3686 


3672 


Vc' 


3638 


3630 


^(S^Si) 


4040 


4072 


r^c(3iSo) 




4043 


V^(43Si) 


4415 


4406 


r/c(4iSo) 




4384 


X2 


3556 


3556 


Xi 


3511 


3505 


Xo 


3415 


3424 


h. 




3516 


X2(2^P2) 




3972 


Xi(23Pi) 




3925 


Xo(23Po) 




3852 


K{2^V,) 




3934 


^S^Dg) 




3806 


^2('D2) 




3800 


^(^Di) 


3770 


3785 


r/c2('D2) 




3799 


^3(2^03) 




4167 


^2(2^02) 




4158 


^(2=^Di) 


4159 


4142 


7702(2^02) 




4158 


X4(^F4) 




4021 


X3(^F3) 




4029 


X2(^F2) 




4029 


r/c3('F3) 




4026 
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TABLE 4. Theoretical and experimental El radiative partial widths of the easily accessible 
1 states, as well as some interesting additional cases (see text). The experimental ^(3770) 
numbers are taken from Ref.|lUj. 



Initial state 


Final state 


^thy. (keV) 


rex.pt. (keV) 


X2 


-fj/ijj 


424 


426 ± 48 


Xi 


'jj/ijj 


320 


288 ±51 


Xo 


7J/^ 


155 


119 ±17 


he 


IVc 


494 




i^' 


1X2 


38 


18.0 ±2.0 




IXi 


54 


23.6 ±2.7 




7Xo 


62 


24.2 ±2.5 


Vc' 


iK 


49 




^(33Si)(4040) 


1X2 


0.5 






1X\ 


0.4 






7Xo 


0.2 






7X2(2='P2) 


14 






7X1(2=^?!) 


39 






7Xo(2='Po) 


54 




7/;(3Di)(3770) 


7X2 


4.9 


< 330 (90% c.l.) 




7X1 


126 


280 ± 100 




7Xo 


405 


320 ± 100 


^(23Di)(4159) 


7X2 


0.8 






7Xi 


14 






7Xo 


27 






7X2(2=^P2) 


5.9 






7Xi(2^Pi) 


168 






7Xo(2=^Po) 


485 






7X2('F2) 


51 
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TABLE 5. Open-charm strong decay modes of the 1 states accessible at BES and CLEO- 
c. A reaction-dependent factor has been removed from the decay subamphtudes in the final 
column of the table, so only the amplitude ratios are physically meaningful. 



State 


Mode 




Fe^pi. (MeV) 


Tthy. (MeV) 


Subamps. 


V^(3770) (3Di) 


DD 






43. 






all 




23.6 ±2.7 


43. 




^(4040) (S^Si) 


DD 

DD* 

D,D, 






0.1 
33. 

8. 






D*D* 






33. 


iPi = +0.056 
^Pi = -0.251 
5Fi = 




all 




52 ±10 


74. 




V^(4159) (2 3Di) 


DD 
DD* 






16. 
0.4 






D*D* 






35. 


ip^ = +0.081 
5p^ = -0.036 
^Fi = -0.141 




D,D, 






8. 






all 




78 ±20 


73. 




^(4415) (4 3Si) 


DD 
DD* 






0.4 
2.3 






D*D* 






16. 


iPi = -0.018 
5p^ = +0.081 
5Fi = 




D,D, 






1.3 






D,D/ 






2.6 






D/D,* 






0.7 


iPi = +0.006 
5p^ = -0.028 
5Fi = 




S+P modes 


m 










all 




52 ±10 
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